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I. INTRODUCTION 



Among the six known quark flavours known to date, the top quark is of particular 
interest: Its large mass implies that very high energies are involved in the production and 
decay of this particle, which in turn allows for tests of the fundamental interactions at these 
high energy scales. Moreover, the interactions of the top quark can be studied in greater 
detail than those of the lighter particles since the top quark essentially behaves like a free, 
but extremely short-lived particle. With a mass of m 175 GeV, the lifetime of the top 
quark is about 5 x 10~^^ seconds. This short lifetime effectively cuts off the long distance 
QCD dynamics. In particular, the top quark polarization is not diluted by hadronization 
and thus becomes an additional observable to test perturbative QCD, or, more generally, 
short distance physics. 

An ideal machine to study the properties of top quarks in detail would be a high- 
luminosity, high-energetic e^e~ linear collider. The physics potential of such a machine is 
described for example in [Q]. We just mention here that at center-of-mass energies in the 
range a/s = 400 — 1000 GeV, an annual yield of the order of 10^ top quark pairs may be 
expected. 

For the process e^e~ — > tiX, the production cross sections for longitudinally and 
transversely 0] polarized top quarks are known to order as- The correlations between 
the spins of top quarks and antiquarks have been studied extensively in leading order [Q. 
The longitudinal spin-spin correlations have also been calculated in next-to-leading order 
(NLO) 0,^. Polarization phenomena in top quark pair production near threshold have 
been investigated in 0. 

A convenient theoretical framework to discuss spin phenomena is the concept of the spin 
density matrix, and the main objective of this paper is to present results for the full spin 
density matrix of the ti system to order a^. This allows for a systematic study of spin effects 
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in e^e^ tiX. For phenomenological applications, our results should be supplemented by 
the decay matrices at NLO for the different t and t decay channels 

An alternative approach to the analysis of spin effects in top quark production and decay 
is the computation of the relevant helicity amplitudes. This was accomplished at next-to- 
leading order in [jlO[, where also a Monte Carlo event generator for the case of semileptonic 
ti decays was constructed. 

The outline of the rest of this paper is as follows. We start in section by introducing 
the spin density matrix formalism and apply it to the reaction e^e~ — ti at leading order. 
In section |Tl|we compute the QCD radiative corrections to the results of section |I[ Section 
[TV] contains numerical results for a variety of spin observables. We exhibit their dependence 
on the cm. energy and on the top quark scattering angle and further study the effects of 
electron beam polarization. 

II. KINEMATICS AND LEADING ORDER RESULTS 

In this section we review some basic kinematics and the concept of the spin density 
matrix formalism. To set up the notation, we start with a closer look at the amplitude for 
the process 

e+(p+)e-(p_) ^ (7*, Z*) t{k,)t{h)X, (II.l) 

where e'(e^) denotes an electron (positron) and t{t) describes a top (anti-) quark with mass 
m. We work in leading order in the electroweak coupling and in next-to-leading order in 
the strong coupling as = g1/{4:Ti). To this order the unspecified rest X can be only a gluon. 
The amplitude for the reaction ( |11.1| ) can be written in the following form: 

(11.2) 
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In ( [II. 21) , s = {p^+p^y, Qt denotes the electric charge of the top quark in units of e = yiyra, 



and gl, gl are the vector- and the axial-vector couplings of a fermion of type /, i.e. 

^,^ = r3^-2g^ sin^^^, and g[ = Ti , (II.3) 

in particular (7^ = — | + 2 8110^ ^w, da = ~| fo^ electron, and gl = \ — Isin^t^vK, gi = ^ 
for a top quark, with -t^w denoting the weak mixing angle. The function x(s) is given by 

1 s 

4 sm vw cos^ Vw s — + zm^i z 

where mz and Tz stand for the mass and the width of the Z boson. (We keep here the 
width of the Z boson because it will be relevant for an application of our results to b quark 
production at the Z resonance.) The amplitudes V^, in ( [II. 2| ) encode the information on 
the decay of the vector boson into the ti and tig final states. In particular they depend on 
the momentum and the polarization of the outgoing particles. Considering only longitudinal 
polarization for the incoming electrons and/or positrons and neglecting the lepton masses 
leads to 



fi\ 



S2 



LPC^^-H^^^ + L^^'^'^i/^J'J (II.5) 
for the square of ( [II. 2| ). The lepton tensors 

The tensors iJ^^p^^^^ describing the decay of a polarized Z boson can be written as 

ttPC(PV) _ „VV ttVV I „AA ttAA , VA+ ttVA+ , VA. ttVA- /tt j\ 

^iJLP — ypc{pv)^iiv ypc{pv)^tiv ~^ ypc{pv)^iJ.y ~^ ypc{pv)^iJ.v ' ^-'--'--'J 

with 

= V,V:, hX = A,A:, and H^^ = V,Al ± A,V: . (II.8) 
The couplings g\ (X G {PC, PV}, Y e {VV, AA, VA+, VA_}) in (^[H^ are given by 



9pc{PV) ~ Qt fpciPV) + 2 (?* R'ex('5) fpc{pv) ~^ 
9pc(PV) = 9^a\x{s)?fpc(PV)'> 

d'pclpv) = ~9aQt^^x{s) fpc'iPV) ~ dlgi \xis)\'^fpc(pv)^ 

9piipv) = i 9iQtlmx{s)fJ,l(Pv)^ (11.9) 



where 



(11.10) 



/II = (1 - \-K){gf + gf) - 2(A_ - X+)glgl, fH = 1 - A_A+, 
= (A_ - X+){gf + gf ) - 2 (1 - \^\^)gl gl fH = A_ - A+, 

PpI = -(1 - >^-K)9l + (A- - K)g:, 

ffy = (1 - A_A+)(?^ - (A_ - A+)^7„^ 

with A- (A-I-) denoting the longitudinal polarization of the electron (positron) beamQ. The 
couplings g^c{PV) formally of higher order in the electroweak couplings. The structure 
"will therefore not be discussed further. For top quark production, where a/s ^ rriz, 
one should set the width Tz of the Z boson to zero for consistency. 

The (unnormalized) spin density matrix for the reaction (11.1) may be defined by 



Paa',/3/3' = '^{t{kt, a)t{kt, a')X\T\e^ {p+, A+)e (p__, A_)) 

{tih, /3')^l^|e+(p+, A+)e-(p_, A_))*, (11.11) 

where a, a', (3, (3' are the spin indices of the outgoing top (anti-) quarks. The sum in ( |11.11| ) 
runs over all unobserved degrees of freedom such as the colour of the outgoing particles or 
the polarization of the emitted gluon. In ( [1I.11| ) one should read the combination aa' {[3(3') 



on the left-hand side as a shorthand notation for a multi- index built from a, a' {(3,(3'). To 
calculate the spin density matrix it is convenient to use a different representation which 
follows immediately from the concept of the density matrix: 



For a right-handed electron (positron), A=p = -|-1. 
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Tr 



1 1 
p- -(]l + St-cr)®-(U + SfO-) 



(11.12) 



Here (sf) is the unit polarization of the top (anti-) quark in the rest frame of the top 
(anti-) quark0, and (Tj are the usual Pauli matrices. With (g) we denote the tensor product 
between the spin space of the quark and the antiquark. Using in (|11.12|) a decomposition of 
the spin density matrix p of the form 



p = a]l(g)]l + B+-cr(g)]l + ]l(g)cr-B + djai (g) a 



31 



(11.13) 



the density matrix can be easily calculated by a comparison of the polarization independent 
parts, terms proportional to (sfj), and terms proportional to SuSij on the left-hand side 
and the right-hand side of (|11.12|) . More precisely we define 



Y,X 



(X G {PC.PV}, Y e {VV,AA,VA+}), with 
1 , . . . 1 



Tr 



■ -(ll + SfCr)®-(U + %.o-) 



1 4 



Y 
fifi 



(11.14) 



(11.15) 



where Nc is the number of colours, and Qx are the couplings as given in ( |11.9| ). For the 
density matrices p§ we use a representation as in (p.l3|) . It is useful to decompose the 
polarizations By' and the spin-spin correlations Cyij further. For the two-parton final 
state it is convenient to write: 

B± = bfp + bfk + bfh, 

Cij = Co5ij + SijkiciPk + cskfc + cghfc) + c^PiPj + c^kikj + CQ{pikj + p^k^) 

+ C7{-pifij + Pjfii) + csikitij + kjUi), (11.16) 



^We define the rest frame of the outgoing top (anti-) quark as the rest system which is obtained 
by a rotation-free Lorentz-boost from the center-of-mass system of the e+e~-pair. 
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with 



P 



|kt| 



n 



p X k 



(11.17) 



IP-I \^t\ |p X k| 

where the three-momenta p and k are defined in e'^e~ cm. system. In ( p.l6| ) we suppress for 
simphcity the additional indices Y,X. For the case of the three-parton final state a similar 
decomposition can be used. A detailed discussion of the properties of p under discrete 
symmetry transformations is given in [|Tl|]. In leading order (0(a^)) the non- vanishing 
entries in the density matrices py read: 



PC 



c 



c 



PC 

oyv 

PC 

4, VV 

PC 

5, VV 



2-/32(1-^2- 



±,PC 
,VA+ 



2prz, 



^PC 
.PC 



-2/32, 



2((l-r)V + /32), <S^ = 2/32z 



cl^y = -2(1 -r)z, 



(11.18) 



±,PV 

vv 



2r, 



~ 2/3^2;, 



PV 

'■VA+ 



4/3z, 



b^^^ = 2{l-r)z, 



c!:(A^=m^-r)z, 



^PV 

-6,1^4 



2/3r, 



where 2; = p-k, /3=yl — Am? / and r = 2m/ ^/s. 

The leading order differential cross section (icro(s(,Sf) is related to the leading order 
density matrix po as follows: 

1 



da{st, St) = — Tr 



Po ■ ^{1 + St ■ a) + St ■ a) 



dRo 



(11.19) 



with 



dR2 = ./.h,.J''!L. {2ny6{p+ +p--kt- k-t 



{2Ttf2k^t (27r)32A;5 
The total cross section for example can be obtained from 



(11.20) 
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2s IGn 



1 , 1 

I dz Tr[po] = -Tra^iVc/? | rfz ((7^^ a^^ + gj,^ a^^), (11.21) 
-1 -1 



yielding the well known result: 

(To = CTptNcP 



^ ^\Vc+P'9^pi)^ with a,, = ^. 



(11.22) 



Within the framework of the spin density matrix formalism it is easy to calculate spin 
observables. For instance, at leading order the polarization of the top quark projected onto 
its momentum direction can be obtained from: 



-1 



Po 



k.f 
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/d2;Tr[po] 
-1 



1 

-1 



^!dz{gl^a^<^y + gj.^a'jS,) 
-1 



VA+ 
PC 



(3-/52)(7]^^ + 2/5V^' 



(11.23) 



where = (g) U is the top quark spin operator. (The spin operator of the top antiquark 
is Sf = 11 ® -y.) As another example consider the following spin-spin correlation, which is 
in leading order proportional to the so-called longitudinal spin-spin correlation studied in 



((k-S,)(k-S,-)) 



1 r f 

J dz Tr ^po ■ (k 



f ®k.f 



JdzTrlpo] 



l {l + (3')gVc+2(3^g^^ 
^{'i-PVpl + 2P^g^^- 



(11.24) 



The examples above show that the spin density matrix formalism enables one to calculate 
efficiently the expectation values of spin observables. A more exhaustive analysis of spin 
observables together with next-to-leading order numerical results will be presented in section 

m 



III. QCD RADIATIVE CORRECTIONS 



The QCD corrections at order to the expectation values of spin observables are given 
by the contributions from one-loop virtual corrections to e^e~ ti and from the real gluon 
emission process e~^e^ —>■ tig at leading order. We first give some details on the computation 
of the virtual corrections. 

Both infrared (IR) and ultraviolet (UV) singularities which appear in the one-loop in- 
tegrals of the virtual corrections are treated within the framework of dimensional regular- 
ization in d = 4 — 2e space-time dimensions. We use the 't Hooft-Veltman prescription 



12 1 to treat the 75 matrix present in the axial vector current part of the vertex correction 



in d dimensions. It is well known that this prescription violates certain Ward identities. 



They are restored by adding a finite counterterm |T3| . The UV singularities are removed by 
appropiate counterterms fixed by on-shell renormalization conditions for the quark. After 
renormalization one obtains UV finite vertex corrections for the vector and the axial vector 
parts of the amplitude to order as- 

The renormalized amplitude still contains an IR singularity which appears as a single 
pole in e and which multiplies - up to a factor - the Born amplitude. This singularity 
is cancelled in infrared safe quantities by a corresponding singularity from the real gluon 
emission process. The latter singularity is obtained from the phase space integration of the 
squared matrix element for e^e^ tig over the region of phase space where the gluon is 
soft. 

The virtual corrections to the density matrix are obtained by first computing the in- 
terference between the renormalized one-loop amplitude and the Born amplitude for given 
polarization vectors s^, Sf and then extracting p™^^'^^ as described in section ^ below equa- 
tion ( [11.12| ). Note that the necessary trace algebra can now be performed in d = 4 dimensions 
without punity. In particular, the projectors (1 + 'j^^t,!)/'^ can be kept in 4 dimensions. 
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We now discuss the contributions from real gluon emission. We isolate the soft gluon 
singularities by splitting the tig phase space into a soft and a hard gluon region. The soft 
gluon region is defined by the condition 

Eg < a;^in^, (III.l) 

where Eg is the gluon energy in the cm. system and Xmin is a sufficiently small quantity. 
The hard gluon region is the complement of the soft region. In the limit where the gluon 
momentum kg goes to zero one can neglect kg in the numerator of Tfi{e~^e~ tig), which 
leads to 



Using ( [ill. 21 ) in the whole soft gluon region leads to the approximation 



/ {S^-^^2E - ^M^-"^- - ttg) ^ Sp,{e^e- ti) ^ (III.3) 



where the soft factor S is given by 



d!^ ^kg y/s , f 2ktkt m 



J 27vY-^2E„ ^ '""^ 2 • 



2 



, . \{hkg){ktkg) {hkgf {ktkg)\ 



27r "r(l -e) \ s J ^ e/3 
2e In(cu) + (1 + (^Li2 (1 - a^) + ^ ln2(cu; 



0(6). (III.4) 



Here, Cp = {N^ - l)/(2A^c), P = y^l - 4mVs, and = (1 - + The scale /i 

is introduced in ( [1II.4| ) to keep the strong coupling constant dimensionless in d dimensions. 



The dependence on fi cancels in the sum of the virtual and soft contributions. 

For finite Xmin, the sum of the contributions from the soft and hard gluon region differs 
from the exact result by terms of order 

^min bcCclllSG of the soft gluon approximation. The 
sum becomes exact for x^iin — 0. With the choice Xmin = 10"^ the systematic error due 
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to this approximation is smaller than one permill in all our numerical results. This can be 
nicely checked by varying x^i^ between, say, 10~^ and 10~^ and numerically extrapolating 
to zero. 

The sum of the virtual and soft contributions to the density matrix p is finite and can 
be written in a compact form as follows: 
We define: 



L = -^C'p^l (2/3 + (1 + f3') Hu)) [in(xLn) - In f ) + 2 



+ (1 + (3^) (4Li2 (1 - c^) + \n\uj) - TT^) |, (III.5) 
and use as further abbreviations 

«: = ^C^, ii = -f^PHu), 4 = (2-/5^)4, 4 = ^4. (III.6) 
Then, 



where the nonvanishing building blocks of po are listed in equation ([II. 181) of section 
II. The matrix p'^^'^^ is also decomposed according to equation (PI.14|) with matrices py'^'^^^ 
expanded like in ( [1I.13| ), ( [II.16| ). The nonvanishing entries of the various matrices py'^^^^ 
that make up p^*^^* read (we suppress here the index "rest" for aesthetic reasons): 
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a^^ = (l + z2Ki, = (1 + = - 2)^3, 

btvv = -^nrPzVT^, 4% = (1 - z'Y^, b^;^f^ = (2(1 + z^) + r(/?2 _ 2)z')i„ 



„PC 
'^4,VV 



= 2(1 + (1 - r)zyi 



„PC 
-■6,VV 



-z{2 - r)h, 



b^;^^ = z{2 - r)i,, = f^^riP^ " ^)VT^, 

clly = -KTTrPVT^, cIXa^ = 2z{r{P^ - 2) + 2)^, 



(IIL8) 



For a given observable, the contributions from gluons with energy Eg > x^i-ay^ /2 are 
calculated by a numerical integration over the hard gluon region of the three-body phase 
space. The spin density matrix p^'^'^'^[e'^e~ tig) for the hard gluon emission process is 
obtained by evaluating the left-hand side of equation {^l. 12|) for X = g. The individual 



matrices py'^'^'^'^ are rather lengthy and we do not list them in this paper. We just mention 
here that instead of the expansion ( [II.16| ) of B^, Cij with respect to p, k, and n that was 
used for the two-parton final state, we found it more convenient for the three-parton final 
state in the hard gluon region to use as basis vectors k(/|k(|, kf/|kf|, and (k^ x kf)/|kt x kf |. 
Note that the matrix p^'^^'^ [e~^ e~ — > ttg) does not contain any singularities and that the 
whole computation can be performed in ci = 4 dimensions. 
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IV. NUMERICAL RESULTS 



In this section we present next-to-leading order results for expectation values of a variety 
of spin observables. For an observable O we use the notation 

„2\ 



(0) = (0)o + ^(0)i + oR 



TT 



etc 



(7 = (To H 0"i + O 

TT 



where a is the total cross section for e^e 



ttX, and 



{O)o = -^ fdR2TT{po-0}, 
an Is J 

dR2 Tr jhm ( 



{Oh 



1 1 

1 1 

(To 2s 



^soft _j_ ^virtual j ^ 



Here, dR2 is given in ( [il.2q ) and 

^ (27r)32A;? (27r)32A;e (27r)32A;0 
We consider the following set of observables: 



{27l)^6{p++P_-kt-kt-kg). 



(iv.i: 



(IV.2) 



(IV3) 





= P ■ 


St, 




= P 


St, 


02 


= k- 


St, 


02 


= k 


St, 


0, 


= n ■ 


St, 




= fi 


St, 


04 


= St 


•St, 








0, 


= P ■ 


(St X St), 










= k- 


(St X Sf), 








Oj 


= n • 


(St X Sf), 










= (P 


■St)(p-Sf) 
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= (k 


SO(k-Sf), 








= (P 


Si)(k-Sf) + (k 


S*)(P 


Sf), 


Oil 


= (P 


Si)(n-Sf) + (n 


S*)(P 


Sf), 


Ol2 


= (k 


SO(n-Si) + (n 


S*)(k 


Sf). 



(IV.4) 



The expectation values (02)0 ^^nd (09)0 are given in analytic form in ( [11.23| ) and ( [11.24| ), 
respectively. 

Several constraints are imposed by discrete symmetries on the expectation values of the 
observables ( [1V.4|) . An unpolarized e^e~ initial state is an eigenstate of the combined charge 
conjugation (C) and parity (P) transformation. CP invariance of the interactions considered 
here then implies {Oi) = (Oi) and {O^) = 0. Further, differences between {O2) and {O2) as 
well as nonzero values for (Oq) and (Oj) can only be generated by the contributions from 
hard gluon emission, since Si^Sf, — kj, and since we have kf = — k^ for a final state 
consisting solely of a ti pair (recall that the three-momenta are defined in the e~^e~ cm. 
system). From invariance under the time reversal operation T it follows that nonzero (O3), 
{O3), {Oq), {Oh), and (O12) can only be generated by absorptive parts in the scattering 
amplitude. To order this means that (Oq) is exactly zero due to CP invariance, while 
{O3) = {O3), (Oil), and {O12) get nonzero, albeit small, contributions from the imaginary 
parts of the one- loop integrals appearing in the virtual corrections (cf. the functions bfyy, 



b^y'X_^_, and c^yy, Cjsva+ equation (|1II.8D ). All the above arguments also hold for the 



case of polarized electrons (and/or positrons), although in that case the initial state has no 
definite CP parity. This is because the net effect of a CP transformation of the initial state 
is Aip — A-t in our formulas, and hence the couplings gy are left unchanged, (cf. (|II.9| ), 

(MM))- 

In Table I we list our results for the expectation values of ( [IV. 4| ) in terms of the quantities 
((9i)o,i as defined in ( |IV.1| ) and ( P^V.2| ). We choose four different cm. energies, namely 
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^/s = 400, 500, 800, and 1000 GeV. The positron beam is always assumed to be unpolarized, 
while for the electron beam the three cases A_ = 0, ±1 are considered. As numerical input we 
use mz = 91.187 GeV, an on-shell top quark mass of m = 175 GeV, and sin^ = 0.2236. 

The table shows that the top quark and antiquark are produced highly polarized and 
also that the spin-spin correlations are large. For example, the polarization^ of the top 
quark projected onto the beam axis at -^i = 500 GeV and for A_ = +1 amounts to 

2(p . St) = 0.8998 - 0.278— = 0.8910, (IV.5) 

71 

where we set = 0.1. As another example, consider the spin-spin correlation (Oio) at 
y/s = 1 TeV, also for A_ = +1: 

((p • St)(k ■ St) + (k ■ St)(p ■ Sf)) = 0.2770 - 0.303— = 0.2674, (IV.6) 

TT 

where we again set = 0.1. 

A global characteristic of all the expectation values of the observables (|IV.4|) is that the 
QCD corrections are quite small. The quantity a^/vr x \{Oi)i/ {Oi)o\ ranges, for nonzero 
((9j)o and (a fixed value of) Og = 0.1 between 1.9 permill (for {O4) at ^/s = 400 GeV and 
all three choices of A_) and 5.3 percent (amusingly also for (O4), but at ^/s = 1000 GeV 
and A_ = -1)5 

To check our calculation, we compared our numerical value for the order correction 
to the total cross section ai with the value one gets by using the analytic formula as 



given for example in and found excellent agreement. Note that the longitudinal spin- 
spin correlation (P^^) studied in p is, at next-to-leading order, not proportional to our 
expectation value (Og): The former would correspond in our notation to the expectation 



The polarization is conventionally defined as two times the expectation value of the spin operator. 
^In leading order, (04) = 1/4, since the reaction proceeds through a single spin-one boson. 
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value 4((kj ■ St)(kf ■ Sf)), which only at leading order is equal to —4(09). To compare our 
results for (P^^), we reproduced Figures 1 and 2 of reference 0| and found agreement. 

We now study the distributions of our expectation values with respect to z, the cosine 
of the top quark scattering angle in the cm. system. These distributions are defined as 
{Oi6{z — z')), i.e. we do not average over z but over all other kinematic variables. 

The distributions (03,ii.i2<^(2; — z')) are not shown, since they can be easily constructed 
from the listed analytic formulas for bf^cj^s of equation ( [111.8| ). We also do not show the 
distribution {O^Slz — z')) , since according to Table I the expectation value (Oj) varies (again 
for a fixed a.^ = 0.1) between the tiny values —0.6 x 10"'^ and —0.3%. 

Figs, la and lb show, to NLO accuracy, the distribution {Oi6{z — z')) at cm. energies 
y/s = 500 GeV and y/s = 1 TeV , respectively, for A_ = 0, ±1. In this and all the following 
plots we set = 0.1. Note that the distribution gets more peaked near z = +1 as the 
cm. energy rises. This feature is less pronounced in the distribution {02S{z — z')) depicted 
in Figs. 2a,b. In Figs. 3 - 6 we show the distributions for different spin-spin correlations, 
namely (04,8,9,io^(-2^ ~ z')). In these figures, the cm. energy is varied between ^/s = 400 
GeV and y/s = 1 TeV, while the electron polarization is set to A_ = 0. The results for other 
choices of A_ do not differ much from the ones shown. This is also reflected in the rather 
weak dependence of the spin-spin correlations (C4,8,9,io) on A_ (cf. Table I). Note that the 
distributions typically rise as z +1. 

To illustrate the impact of the 0{as) corrections, we plot in Figs. 7-10 the "i^-factors" 

^•<^' = {0.5(z-z%, (■^■'> 

for i = 1 (Figs. 7a,b), and i = 4, 8, 9 (Figs. 8,9,10). The i^-factors show a strong dependence 
both on the cosine of the scattering angle and on the cm. energy. They vary between 0.88 
and 1.04. 
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V. CONCLUSIONS 



The production of top quark pairs in e'^e~ annihilation involves a variety of spin phe- 
nomena. We have performed a systematic study of these effects to order as and including 
beam polarization effects using the spin density matrix formalism. Apart from a signifi- 
cant polarization of the top quarks and antiquarks, the spins of t and t are also strongly 
correlated. The QCD corrections to the leading order results for the expectation values 
of all spin observables considered are at the percent level or smaller. The spin effects in 
the tt production will manifest themselves in the angular distributions of the t and t decay 
products. For a phenomenological analysis of these angular distributions, one can combine 
the results presented in this paper with spin decay matrices computed to next-to-leading 
order accuracy for the different t and t decay modes. 
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TABLE CAPTION 



Table I. Expectation values of the observables listed in (|IV.4| ) in terms of the quantities 
(Ci)o,i as defined in ([IV. 1| ) and ( [IV. 2|) for different cm. energies, A+ = and A_ = 0,±1. 
For the expectation values not listed in the table we have, as discussed in the text, ((5i,3) = 
(Oi,3), and (Os.e) = 0. 
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FIGURE CAPTIONS 



Fig. 1. Expectation value ( Oi6{z — z') ) to order for a fixed value = 0.1 and A+ = 0. 
In la (lb) the cm. energy is set to ^/s = 500 GeV {^/s = 1 TeV). The solid line is the 
result for A_ = 0, the dashed line for A_ = —1, and the dotted hne for A_ = +1. 
Fig. 2. Same as Fig. 1, but for ( 025{z - z')). 

Fig. 3. Expectation value ( O^S^z—z') ) to order for a fixed value = 0.1, A+ = A_ = 0, 
and cm. energies a/s = 400 GeV (dashed line), y/s = 500 GeV (solid line), y/s = 800 GeV 
(dotted line), and y/s = 1000 GeV (dash-dotted line). 
Fig. 4. Same as Fig. 3, but for ( OsS{z - z')). 
Fig. 5. Same as Fig. 3, but for ( 0(j5{z — z')). 
Fig. 6. Same as Fig. 3, but for ( Oiq5{z — z')). 

Fig. 7. The function Ki{z) defined in equation ( |1V.7D . In 7a (7b) the cm. energy is set 
to y/s = 500 GeV {^/s = 1 TeV). The solid line is the result for A_ = 0, the dashed line for 
A- = —1, and the dotted line for A_ = +1. 

Fig. 8. The function K^lz) defined in equation ([IV. 7] ) for A^ = A_ = and cm. energies 
v^i = 400 GeV (dashed line), = 500 GeV (sohd line), ^/s = 800 GeV (dotted line), and 

= 1000 GeV (dash-dotted line). 
Fig. 9. Same as Fig. 8, but for Ks{z). 
Fig. 10. Same as Fig. 8, but for Kg{z). 
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